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INTRODUCTION 


The  NAVSTAR  Geodetic  Receiver  System  (NGR S)  records  both  pseudorange  and 
Doppler  observations  from  the  Global  Positioning  System  (GPS)  satellites. 
Observations  of  pseudorange  are  available  each  6  sec  and  a  Doppler  count  each 
minute.  In  order  to  use  these  observations  to  do  precise  point  positioning 
with  the  receiver  or  to  do  satellite  orbit  calculations,  the  time  of  each  ob¬ 
servation  must  be  known  with  high  accuracy  in  the  GPS  time  system.  At  sites 
remote  from  the  GPS  Master  Control  Station  (MCS) ,  GPS  time  must  be  obtained 
indirectly,  since  the  satellites  do  not  transmit  time  signals  that  are  synchro¬ 
nized  with  those  kept  at  the  MCS.  The  satellites  do,  however,  broadcast  cor¬ 
rection  coefficients  in  the  navigation  message.  The  NGRS  method  of  converting 
local  time  to  GPS  time  is  described  in  the  following  sections. 


USE  OF  TIME-TAGS 


The  Doppler  count  start  and  stop  times  are  controlled  by  the  local  clock 
minute  epochs  and  are  therefore  time-tagged  with  these  epochs.  These  epochs 
are  not,  however,  synchronized  with  the  GPS  time  system  and  must  be  corrected 
in  some  manner  so  that  they  also  reflect  true  GPS  time.  The  accuracy  to  which 
this  correction  must  be  made  can  be  estimated  in  the  following  way. 


A  worst-case  relative  velocity  between  a  satellite  at  GPS  altitude  and  a 
station  on  the  earth  is  about  1400  m/sec. *  If  1  cm  is  taken  as  an  allowable 
range  error,  then  the  time  required  for  the  satellite-ground  station  range  to 
change  by  1  cm  is 


At 


0.01 


7.1  ysec 


Therefore,  if  time  tags  on  the  NGRS  observations  are  corrected  to  GPS  time  with 
an  accuracy  of  about  ±7  psec,  the  error  in  the  calculated  satellite  range  will 
be  less  than  1  cm. 


The  satellite  position  is  obtained  from  a  post-fit  ephemeris  that  lists 
positions  at  specific  GPS  times.  Positions  intermediate  to  those  given  are 
obtained  by  interpolation.  Thus,  it  is  clear  that  the  correct  GPS  time  is 
essential,  if  the  satellite  position  is  to  be  calculated  accurately  at  the  time 
of  each  observation.  Once  calculated,  it  is  used  in  one  of  two  ways  depending 
upon  the  problem  being  studied. 

If  improvement  in  the  satellite  orbit  is  the  object,  then  the  calculated 
satellite  position  is  used  as  the  nominal  to  form  differences  with  the  observa¬ 
tions.  It  is  important  that  this  nominal  be  time-tagged  correctly  from  obser¬ 
vation  to  observation  so  that  a  minimum  amount  of  "noise"  from  the  calculated 
satellite  position  is  introduced  into  the  computation.  If  absolute  positioning 
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of  the  receiver  is  the  yoa t ,  i_a^n  the  calculated  satellite  positions  are  treat¬ 
ed  as  the  truth,  and  the  site  posi’-ion  is  computed  from  them  by  using  the  ob¬ 
servations.  Clearly  an  error  in  satellite  position  will  translate  into  an 
error  in  the  computed  receiver  position. 


METHOD  OF  TIME-TAG  CORRECTION 


The  satellites  transmit  time  tags  as  part  of  the  GPS  navigation  message.^ 
These  time  tags  are  keyed  to  agree  with  the  pseudorandom  code  sequence  that  is 
continuously  transmitted.  However,  these  time  tags  are  not  aligned  with  GPS 
time,  and,  in  fact,  each  satellite  clock  drifts  by  different  and  time  variable 
amounts.  This  is  a  consequence  of  the  fact  that  each  satellite's  frequency 
reference  is  necessarily  allowed  to  run  free.  This  situation  is  tolerable 
because  the  time  offset  between  each  satellite  and  GPS  time  is  predicted  in 
advance,  and  corrections  to  the  transmitted  time  can  be  computed  by  the  indivi¬ 
dual  user,  if  he  is  able  to  decode  the  satellite  navigation  message.  The  cor¬ 
rections  to  the  apparent  satellite  time  are  designated  ts  in  Figure  1.  The 
other  offsets  that  are  relevant  to  this  problem  are  also  defined  in  that 
figure. 


- ► 

GPS  TIME 


GPS  LOCAL  SATELLITE  LOCAL  TIME 

EPOCH  EPOCH  EPOCH  AT  WHICH 

SATELLITE  EPOCH  RECEIVED 


rr  =  PROPAGATION  DELAY 
o  EPOCH  OFFSET  (  LOCAL  SATELLITE  I 
•„  -  OBSERVATION  OF  PSEUDO  RANGE  +  ,r 
r(  =  EPOCH  OFFSET  (GPS -LOCAL) 

't  EPOCH  OFFSET  (  GPS  SATELLITE)  =  ra  *  rf 
A  VALUE  FOR  r,  IS  FOUND  IN  THE  SATELLITE  DATA  MESSAGE 


Figure  1.  Diagram  Illustrating  the  Relation  Between  Epochs 
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The  significance  of  the  events  illustrated  in  Figure  1  is  the  following. 

At  a  particular  satellite  epoch,  the  signal  indicating  that  particular  epoch  is 
transmitted.  This  signal  is  received  xr  sec  later  by  the  NGRS .  In  order  to 
calculate  the  geometric  range  corresponding  to  xr ,  the  satellite  position  at 
transmission  and  the  receiver  position  at  reception  are  both  needed  in  a  common 
time  system.  The  satellite  position  at  transmission  can  be  obtained,  if  the 
satellite  epoch  can  be  converted  to  the  corresponding  GPS  time.  The  NGRS  posi¬ 
tion  at  reception  can  be  obtained,  if  the  local  epochs  can  be  converted  to  the 
corresponding  GPS  time.  As  stated  previously,  the  correction,  xs,  to  the 
satellite  epoch  is  available  from  the  individual  satellites'  navigation  mes¬ 
sage.  The  correction  to  the  local  clock,  x^,  involves  some  additional  com¬ 
putations  using  the  observed  pseudorange  xp. 

Using  the  notation  from  Figure  1,  the  GPS  time  of  transmission  is  tgpg  + 
xs,  and  the  locally  generated  epoch  that  corresponds  to  the  satellite  epoch 
is  t^pg  +  tj.  The  GPS  time  of  reception  of  the  transmitted  epoch  is 
CGPS  +  xl  +  Tp*  *n  or<^er  to  evaluate  x^,  we  can  write  an  equation  that 
states  that  the  time  of  transmission  plus  the  geometric  progagation  delay,  xr, 
equals  the  time  of  reception. 

fcGPS  +  Ts  +  Tr  =  hsPS  +  TZ  +  Tp  ( 1  * 

The  local  offset  must  then  Tj  =  is  t  Tr  -  ip,  where  ts  comes  from  the 
navigation  message,  xr  is  calculated  from  the  ephemeris,  and  xp  is  the  pseudo¬ 
range  observation  expressed  in  seconds. 

Calculation  of  rr  requires  the  satellite  position  at  the  time  of  trans¬ 
mission  tGps  +  ts  and  the  receiver  position  at  the  time  of  reception  tGpg  +  is 

+  xr •  However,  since  ir  is  not  yet  known,  we  can  take  the  time  of  reception 
equal  to  the  time  transmission  and  calculate  an  approximate  value  for  Tr .  Then 
we  car.  use  this  xr  to  refine  the  approximation  indefinitely  or  until  the  change 
in  xr  after  each  iteration  becomes  insignificant. 

All  Doppler  count  intervals  recorded  by  the  NGRS  are  adjusted  to  begin  and 
end  at  local  1-min  epochs.  Since  the  local  clock,  like  the  individual  satel¬ 
lite  clocks,  is  free  running,  these  epochs  need  to  be  converted  to  GPS  time. 
Since  we  have  calculated  the  local  offset  i£  from  the  pseudorange  data,  it  is 
a  straightforward  process  to  make  this  correction.  Adopting  the  convention  in 
Figure  1,  the  local  epoch  occurs  after  the  GPS  epoch.  Therefore,  the  GPS  time 
of  the  local  epochs  is  the  indicated  local  time  plus  x^. 

‘xSPS  =  fcLoca  1  +  Tl  <2) 

We  assume  in  doing  this  that  the  indicated  local  epoch  and  the  GPS  epoch  are 
identical.  All  Doppler  time  tags  should  be  corrected  by  the  value  of  x^ 
current  at  the  epoch  of  interest  because  x^  itself  is  a  slowly  varying  func¬ 
tion  of  time. 
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AN  EXAMPLE  OF  LOCAL  OFFSET 


Some  pseudorange  observations  from  the  NGRS  were  made  available  for  a  1 5- 
day  span,  in  order  to  test  the  formulation  described  in  the  previous  section. 
The  reference  ephemeris  for  this  span  was  obtained  from  the  least-square  batch 
fit  to  monitor  station  data  that  is  performed  routinely  at  the  Naval  Surface 
Weapons  Center  (NSWC)  for  the  GPS  Master  Control  Station. 3  The  satellite  clock 
corrections  were  obtained  directly  from  the  satellite  navigation  message. 

Figure  2  is  a  plot  of  the  computed  offset  -t^.  Data  from  four  satellites  are 
plotted  together  in  this  figure  and  illustrate  that  the  computed  local  offset 
is  independent  of  satellite  to  within  a  tolerance  of  about  ±100  nsec.  Since 
the  pseudorange  observations  used  in  this  computation  were  single  frequency 
(Lj),  the  systematic  spread  seen  in  the  figure  is  probably  primarily  due  to  the 
variation  of  ionospheric  delay  with  elevation  angle.  If  the  observations  with 
elevation  angles  less  than  45°  are  deleted  from  the  plot,  as  in  Figure  3,  the 
spread  is  considerably  reduced. 

At  NSWC,  SV-8  passes  within  10°  of  the  zenith  once  a  day.  Data  from  this 
satellite  are  plotted  in  Figure  4  when  the  elevation  angle  is  75°  or  greater. 

At  this  high  elevation,  the  ionospheric  contribution  to  the  L^  range  would  be  a 
minimum,  and  the  observations  would  be  collected  at  nearly  the  same  time  each 
day.  At  this  time  of  the  year,  this  high  elevation  pass  occurred  at  about 
4.4  hr  GPS  or  23.4  hr  EST.  Table  1  lists  the  coefficients  of  a  linear  fit  to 
this  data  and  also  to  all  the  data  plotted  in  Figure  2.  The  root  mean  square 
(RMS)  residual  of  a  linear  fit  to  the  high  elevation  data  over  the  observation 
span  was  11.5  nsec,  which  is  approximately  one  quarter  the  RMS  of  all  the  data 
over  the  entire  span. 

In  effect,  a  time  transfer  has  been  accomplished  with  quite  high  accuracy 
over  more  than  10  days.  The  time  drift  between  GPS  and  the  local  reference  can 
also  be  determined  from  Table  1.  The  fit  parameter,  a^  indicates  the  drift  to 
be  about  66  nsec  per  day  or  0.000764  nsec  per  second.  The  quadratic  term  seems 
to  be  quite  small  and  consequently  has  not  been  evaluated. 

The  deviations  of  the  observations  from  the  linear  fit  are  plotted  in 
Figure  5  versus  the  elevation  angle  of  the  observation.  The  fit  parameters 
used  to  form  these  residuals  were  taken  from  the  fit  to  the  high-elevation  SV-8 
data  only.  The  propagation  path  length  would  be  expected  to  increase  at  lower 
elevations  due  to  the  effects  of  the  ionosphere;  consequently,  the  differences 
of  all  data  with  the  SV-8  fit  should  show  positive  residuals.  The  zero  residu¬ 
al  line  in  Figure  5  should,  therefore,  be  an  approximate  lower  bound  to  the 
offset  residuals.  Figure  b  shows  that  this  is  not  tne  case.  Figure  b  shows 
that  there  are  offset  residuals  both  above  and  below  the  zero  residual  line. 

The  observed  scatter  must  be  due  to  variations  in  the  reference  ephemeris.  An 
error  source  such  as  this  would  produce  residuals  that  were  independent  of 
elevation  angle.  There  does  exist  a  preference  for  positive  residuals  at  low 
elevations.  This  is  probably  due  to  the  expected  propagation  effects  and  also 
to  the  inclusion  of  data  from  satellites  not  observed  at  nigh  elevations. 
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Figure  3.  Clock  Offset,  NGRS-GPS  (nsec)  vs  Date,  SV-5,  6,  8,  9;  NSWC 
Post-fit  Ephemeris,  Elevation  >45°,  NSWC  Dahlgren,  Days  182-197  1980 
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4.  Clock  Offset,  NGRS-GPS  (nsec)  vs  Date,  NSWC  Post-fit  Ephemer 
Elevation  >75°,  GPS  SV-8;  NSWC  Dahlgren,  Days  182-197  1980 
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Figure  5.  Clock  Offset  Residual  (nsec)  vs  Elevation  (°) , 
NSWC  Dahlgren,  Days  182-197  1980 
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Table  1.  Linear  Fit  to  Local-GPS  Clock  Offset 


Ephemer is 

*0 

(nsec) 

al 

(nsec/day) 

RMS 

(nsec) 

PTS 

NSWC 

e  >  75° 

-9438.1 

66.0486 

11  .47 

16 

all  data 

-9237.3 

65.0626 

47.82 

414 

Broadcast 

e  >  75° 

-9543.4 

66.5893 

10.32 

16 

all  data 

-9221  .0 

65.0882 

37.12 

376 

Figure  5  illustrates  that  the  reference  trajectories  of  the  four  satellites 
observed  are  consistent  to  within  about  ±60  nsec  (±18  m) ,  if  the  low  elevation 
bias  toward  positive  residuals  is  neglected.  The  consistency  is  likely  to  be 
better  than  ±60  nsec  because  these  residuals  include  the  predicted  satellite 
clock  correction.  This  clock  correction  is  derived  independently  of  the  NSWC 
satellite  ephemeris.  The  two,  therefore,  cannot  be  expected  to  be  as 
consistent  as  they  might  be,  if  both  clock  correction  and  ephemeris  were 
obtained  simultaneously  from  the  same  data. 


The  time  drift  described  above  can  also  be  interpreted  as  a  fractional 
frequency  deviation  In  order  to  show  this,  the  frequency  standard  can  be 

represented  by  two  equivalent  expressions: 


(t)  =  4>0  + 

2it  fx  (t) 

(3) 

*2  (t*  =  ‘t’o  + 

2tiC  (f)t 

(4) 

In  Equation 
is  the  time 


(3),  the  time  argument,  t,  is  a  function  of  the  parameter,  t,  which 
argument  of  Equation  (4).  For  our  purposes  t ( t )  will  be  defined  as 


t (t)  =  aQ  +  (1+aj )  t 

where  a0  is  the  time  offset  at  t  =  o  and  at  is  the  time  drift  (sec  per  sec). 
This  form  has  the  following  properties: 


t (0)  =  aQ  and  ±(t)=a0+t  ifa1=o 
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Under  the  Latter  condition,  •  n.n-ks  t  exactly  and  differs  only  by  aQ. 

In  t  ,  (  t ) ,  the  frequency  £  is  a  function  of  the  frequency  used  in  <J>j(t), 
Let  £  be  ottset  sliahtly  from  t . 

£  ( t )  =  t  +  At 


Since  these  two  $  expressions  are  equivalent",  they  <-an  be  equated  and  solved 
for  At. 

t  )  -  iji  t  r  i 


/!nt[a,,  +  U  +a 


ii  i  t  *  t  » 


Af  o 

■r  =  —  +  rti 

The  term  -pi  becomes  sm.i  i  i  if  ter  i  no  ’  1  .t>-  ,  so  if  is  [>ermissible  to  take  the 

offset  to  be  zero,  or  r  larqe.  nim  -himi  ftins,  the  fractional 

frequency  deviation  4^-  is  »«qu»v  i  !•••>«  *  *  i  me  lr  1  t  r  . 

a  |  -  4-^-  : , o4  x  y'-*  * 

This  number  represents  the  fractional  frequency  deviation  between  the  local 
cesium  reference  and  the  ensemble  tnat  defines  the  GPS  time  system.  The  speci¬ 
fication  for  the  ionq-term  fractional  frequency  deviation  of  the  local  cesium 
is  ±7  x  10-1-.  This  observed  value  is  well  within  the  specification. 


A  SECOND  EXAMPLE 


The  same  data  span  was  reevaluated  for  the  local  offset  usinq  the  same 
formulation  as  previously  described.  In  this  case,  however,  the  satellite 
trajectory  was  computed  usinq  the  broadcast  ephemeris  instead  of  a  post-fit 
ephemens.  This  technique  allows  more  timely  results  to  be  obtained  because 
all  the  information  needed  is  contained  in  the  data  as  recorded  by  NGRS  from 
the  satellite  NAVUATA  messaqe. 


Eiqure  b  is  the  plot  of  the  computed  offset,  —  x ^ ,  from  all  data.  It  shows 
a  few  more  wild  points  than  the  cor respondinq  fiqure  from  the  NSWC  post-fit 
ephemeris,  but  this  may  be  due  to  our  implementation  of  the  alqorithm  and  not 
intrinsic  to  the  broadcast  ephemeris  itself.  Otherwise,  the  aqreement  between 
the  two  ephemendes  is  excellent.  Eiqure  7  is  a  plot  usinq  the  data  from  SV-8, 
which  was  from  the  portion  of  the  pass  whose  eLevation  anqle  was  qreater  than 
75°.  Table  1  indicates  that  the  solution  to  the  two  sets  of  data  differed  by 
only  5.3  nsec  in  aQ  and  oy  0.5  nsec  per  day  in  a^.  The  RMS  computed  from  all 
data  was  10  nsec  less  for  the  broadcast  ephemeris  when  compared  to  the  NSWC 
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ephemeris,  but  this  was  after  all  observations  yreater  than  three  times  the  kM., 
were  eliminated  from  both  sets  of  data.  No  observations  were  eliminated  from 
the  SV-8  fits  where  the  elevation  angle  was  yreater  than  7b°. 

When  the  fit  to  the  high-elevation  data  is  used  to  form  residuals  with  the 
rest  of  the  observations  evaluated  using  the  broadcast  ephemeris,  the  result  is 
quite  different  from  that  obtained  using  the  NSWC  ephemeris.  In  this  case, 
illustrated  in  Figure  8,  bias  toward  positive  residuals  at  lower  elevation 
angles  is  quite  apparent.  This  would  seem  to  imply  that  the  random  component 
due  to  ephemeris  error  is  somewhat  less  for  the  broadcast  ephemeris.  If  the 
plot  is  used  as  a  guide,  an  estimate  of  the  scatter  is  about  ±40  nsec  (±12  m). 
Since  the  same  satellite  clock  corrections  are  used  in  both  cases  discussed,  a 
comparison  of  Figures  5  and  8  is  an  indication  of  the  relative  consistency  of 
those  clock  corrections  with  each  of  the  ephemerides.  As  expected,  the  broad¬ 
cast  ephemeris  fits  the  broadcast  clock  corrections  better  than  the  NSWC 
ephemeris  does.  This  result  would  follow  from  the  methods  of  computation  em¬ 
ployed  in  each  case.3 

CONCLUSIONS 


The  comparison  of  the  broadcast  ephemeris  with  the  NSWC  post-fit  ephemeris 
indicates  that  the  algorithm  to  implement  the  broadcast  ephemeris  is  being 
applied  correctly.  The  use  of  the  broadcast  ephemeris  is  an  advantage  because 
it  can  be  recorded  on-line  at  the  same  time  the  observations  are  made.  As  a 
by-product  of  this  checkout  procedure,  an  indication  of  the  relative  scatter  of 
local  clock  residuals  computed  by  using  the  two  satellite  ephemerides,  but  the 
same  observations  and  satellite  clock  corrections,  have  been  made.  These 
residuals  indicate  that  the  broadcast  ephemeris  performance  is  slightly  better 
than  the  post-fit  ephemeris  when  used  in  conjunction  with  the  broadcast  clock 
corrections. 

Of  primary  significance  is  the  behavior  of  the  receiver,  peripherals,  and 
local  frequency  standard.  No  significant  problems  were  experienced  throughout 
this  14-day  period,  and  the  data  received  were  of  high  quality.  An  indication 
of  the  data  quality  is  illustrated  by  the  plots  in  Appendix  A. 
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Figure  7.  Clock  Offset,  NGRS-GPS  (nsec)  vs  Date,  GPS  Broadcast  Ephemeris, 
Elevation  >75°,  GPS  SV-8,  NSWC  Dahlgren,  Days  182-797  1980 
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Figure  8.  Clock  Offset  Residual  (nsec)  vs  Elevation 
Dahlgren,  Days  182-197  1980 
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APPENDIX  A 


The  data  contained  in  this  appendix  was  compiled  to  perform  diagnostic 
checks  on  the  operation  of  the  NGRS  hardware.  During  the  182-197-day  span 
considered  in  this  report,  performance  of  the  system  was  quite  satisfactory  as 
the  plots  that  follow  illustrate. 

In  the  NGRS,  range  and  Doppler  observations  are  obtained  independently, 
pseudorange  data  from  L ^  only,  and  Doppler  from  L ^  and  L2  simultaneously.  The 
pseudorange  data  rate  was  one  observation  each  6  sec;  the  Doppler  rate  was  one 
observation  each  60  sec.  In  order  to  produce  the  plots,  data  from  10  to  20  nun 
of  continuous  observations  were  smoothed  by  fitting  the  pseudorange  with  a 
fourth-degree  polynomial,  and  the  range  differences  with  a  third-degree  polyno¬ 
mial.  The  RMS  of  the  residuals  between  the  fit  and  the  observations  is  usea  as 
the  primary  indicator  of  performance. 

Fourth  degree  fits  to  pseudorange  and  third-degree  tits  to  range  differ¬ 
ence  over  20  min  produce  residuals  that  have  tew  systematic  trends  remaining. 
The  residuals  are  indicative  of  the  observational  noise  from  tne  various  con¬ 
tributing  sources  such  as  preamplifier  noise,  receiver  phase  jitter,  frequency 
standard  fluctuations,  and  measurement  error  sources.  Results  obtained  from 


examination  of  the  plots 

are 

summarized  in 

Table  A- 1 . 

Table 

A—  1  . 

Average  RMS 

of  Residuals  (cm) 

SV:  5 

6 

8 

9 

Pseudorange 

RMS 

<v 

(L1  ) 

(v 

76.9 

66.1 

68.8 

76.6 

Range  Difference 

RMS 

3.5 

3.2 

4.1 

6.3 

Range  Difference 

RMS 

3.5 

3.2 

4.1 

6.3 

Range  Difference 

RMS 

(corrected ) 

3.7 

3.4 

4.2 

6.4 

Range  Difference 

RMS 

( ionosphere ) 

0.8 

0.9 

0.7 

0.9 

It  is  interesting  to  note  that  SV-9,  whose  reference  frequency  standard  is 
a  cesium,  indicates  a  higher  RMS  residual  than  the  other  three  satellites  whose 
frequency  standards  are  rubidiums.  The  NGRS  reference  is  an  HP-5061-004  high- 
performance  cesium.  The  3-cm  residual  RMS  is  probably  due  to  this  local  refer¬ 
ence.  Experiments  with  other  cesiums,  rubidiums,  and  quartz  oscillators  have 
caused  this  floor  RMS  to  change  as  would  be  expected  for  the  other  types  of 
standards.  Note  that  the  ionospheric  correction  RMS  is  below  1  cm.  In  calcu¬ 
lating  this  number,  much  of  the  correlated  noise  contribution  contained  in  Lj 
and  L2  is  removed.  For  this  reason  1  cm  is  more  representative  of  the  system 
noise  from  uncorrelated  sources. 

A  set  of  12  plots  is  used  to  judge  the  performance  of  the  system  over  the 
time  span  of  interest.  A  set  is  generated  for  each  sateLLite  observed  so  that 
satellite-related  anomalies  are  readily  apparent. 
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1.  20— MIN  PSKHDORANGK  K MS  rv  VS  PATH 

Toe  RMS  of  the  residuals  ot  a  f our th-deqree  polynomial  fit  to  2o-mn 
BP(jmen:s  ot  b-sco  psenrioran.j  es  are  displayed  in  the  first  plot.  As  was 
noted  in  Table  A-1 ,  the  RMS  from  all  four  satellites  <jive  about  the  same 
results,  and  there  is  little  variation  with  date.  The  7<j-  to  80-cm  RMS 
represents  a  2.7-nsec  jitter  in  the  measurement.  The  vertical  spread  in  the 
RMS  is  due  to  some  extent  to  variations  in  the  elevation  a no  Le  as  is 
illustrated  in  the  next  plot. 
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Figure  A-1.  20-i-lin  Pseudoranye  RMS  (cm)  vs  Hite,  GPS  SV-S, 
tJSWC  Uaiilyren,  Days  182-1 D7  1‘JUU 
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Figure  A-3.  20-Mi.n  Pseudorange  RMS  (cm)  vs  Date,  GPS  SV-8, 

NSWC  Dahlgren,  Days  182-197  1980 
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2.  20-MIN  PSEU1XJRANGK  RMS  (cm)  VS  ELEVATION  lJJ 

In  this  plot,  the  pseudorange  residuals  are  plotted  versus  elevation 
angle.  They  show  a  sliqht  tendency  to  increase  at  lower  elevations.  Tins  is 
best  illustrated  by  SV-8  because  data  are  available  trom  20  to  80°  elevations. 
This  effect  could  have  two  causes.  Variations  in  the  propagation  path  due  to 
ionospheric  or  tropospheric  irregular ities  would  be  more  apparent  at  lower 
elevation  angles;  also  the  signal  strenqtn  will  decrease  sliqhtly  due  to  in¬ 
creased  ranqe  and  reduced  antenna  qain. 
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Figure  A-5.  20-Min  Pseudorange  RMS  (cm)  vs  Elevation  (°), 
GPS  SV-5 ,  NSWC  Dahlgren ,  bays  182-197  1980 
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Figure  A-6 .  20-Min  Pseudorange  RMS  (cm)  vs  Elevation  (*) 
GPS  SV-6 ,  NSWC  Dahlgren,  Days  182-197  1980 
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Figure  A-8.  20-Min  Pseudorange  RMS  (cm)  vs  Elevation  (°) , 
GPS  SV-9,  NSWC  Dahlgren,  Days  182-197  1980 
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RMS  (cm)  VS  DOPPl.RK  (Hz) 


.  20-MIN  Lj  RANGH  Lj I F 1-  i.e  ......  i, 

The  RMS  of  the  residuals  ot  a  third-decree  polynomial  fit  to  20-mu:  neq- 
ments  of  bO-sec  L ^  ranee  differences  (derived  from  Doppler)  is  displayed  versus 
the  Doppler  frequency.  A  positive  Doppler  frequency  indicates  approach.  There 
is  little  variation  with  Doppler  frequency,  and  all  residuals  are  clustered 

around  3.5  cm  except  for  SV-9.  This  can  lie  explained  by  the  knowledqe  that  the 

frequency  standard  beinq  used  in  SV-9  was  a  cesium  not  a  rubidium  as  in  the 

other  satellites.  Since  the  other  three  satellites  show  little  difference,  the 

apparent  3.5-cm  floor  is  probably  due  to  the  local  frequency  standard. 
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Fiqure  A-9.  20-Min  Lj  Ranqe  Difference  RMS  (cm)  vs  Doppler  (!!:•.), 
GPS  SV-5 ,  NSWC  Uahlqren,  Days  182-197  198U 
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Figure  A- 1 2 .  20-Min  Lj  Range  Difference  RMS  (cm)  vs  Doppler  (Hz) 
GPS  SV-9,  NSWC  Dahlgren,  Days  182-197  1980 
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4.  20-MIN  L,  RANGE  DIFFERENCE  RMS  (cm)  VS  ELEVATION  (°)  1 

1  .  ] 

The  same  range  difference  residuals  appearing  in  plot  3  are  presented 

versus  elevation  angle  in  plot  4.  The  data  from  SV-8  indicates  that  there  is 

little  or  no  variation  from  20  to  80°.  Only  SV-6  suggests  that  low-elevation  j 

angles  and  larger  RMS  residuals  correlate.  j 
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Figure  A-16.  20-Min  Lj  Range  Difference  RMS  (cm)  vs  Elevation  (°), 
GPS  SV-9 ,  NSWC  Dahlgren,  Days  182-197  1980 
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Figure  19.  20-Min  L2  Range  Difference  RMS  (cm)  vs  Doppler  (Hz) 
GPS  SV-8 ,  NSWC  Dahlgren,  Days  182-197  1980 
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Figure  A-20.  20-Min  L ^  Range  Difference  RMS  (cm)  vs  Doppler  (Hz) 
GPS  SV-9 ,  NSWC  Dahlgren,  Days  182-197  1980 
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20-MIN  CORRECTED  RANGE  DIFFERENCE  RMS  (cm)  VS  DORRLER  (Hz) 


The  RMS  of  the  ionospheric  corrected  range  difference  residuals  shows  that 
combining  the  Lj  and  L0  observations  to  obtain  the  correction  does  not  degrade 
the  observational  accuracy  to  any  great  extent.  This  result  is  in  contrast  to 
the  corrected  pseudorange  observation  error. 
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Figure  A-21 .  20-Min  Corrected  Range  Difference  RMS  (cm)  vs  Doppler  (Hz), 
GPS  SV-5 ,  NSWC  Dahlgren,  Days  182-197  1980 
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7. 


20-MIN  CORRECTED  RANGE  DIFFERENCE  RMS  (cm)  VS  DATE 


The  RMS  of  the  ionospheric  corrected  range  difference  residuals  is  also 
plotted  versus  date,  in  order  to  identify  any  time  variations  that  may  have 
occurred.  Satellites  5,  6,  and  8  show  little  change;  however,  SV-9  may  have 
improved  slightly  after  day  190. 
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Figure  A-25.  20-Min  Corrected  Range  Difference  RMS  (cm)  vs  Date, 
GPS  SV- 5,  NSWC  Dahlgren ,  Days  182-197  1980 
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8.  20-MIN  IONOSPER1C  CORRECTION  TO  L,  RMS  (cm)  VS  DOPPLER  (Hz) 

The  RMS  of  the  ionospheric  corrections  is  plotted  versus  the  Doppler 
frequency  in  plot  8.  Satellites  5,  8,  and  9  are  well-behaved  with  RMS  values 
near  1  cm.  The  SV-5  and  SV-8  are  particularly  impressive  with  most  values  less 
than  1  c.u.  The  SV-9  and  SV-6  demonstrate  a  wider  scatter  whose  origin  must  be 
due  to  a  pecularity  of  these  satellites  or  propagation  paths.  Since  the 
satellite  observations  were  sequential  and  interleaved,  the  receiver  cannot  be 
the  cause  of  the  observed  differences. 
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Figure  A-29.  20-Min  Corrected  Range  Difference  RMS  (cm)  vs  Date, 
GPS  SV-5,  NSWC  Uahlgren,  Days  182-197  1980 


'i.  2U-MIN  LONOSf KKRK-  CORRECTION  TU  Lj  RMS  (cm)  VS  DATE 

The  same  residuals  from  plot  8  are  plotted  versus  date  in  plot  9.  There  i 
no  particular  variation  apparent  with  time  and  so  the  scatter  cannot  be  due  to 
a  temporal  change  in  the  hardware. 
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Figure  A-33.  20-Min  Ionospheric  Correction  to  Lj  RMS  (cm)  vs  Cate, 
GPS  SV-5 ,  NSWC  Dahlgren,  Days  182-197  1980 
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Kiijuro  A-3b.  ^O-Min  Ionospheric  Correction  to  1,^  KMS  (cm)  vs  Date 
GKS  SV-H,  NSWC  bahl<jren,  bays  1  37  1360 
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Figure  A-36.  20-Min  Ionospheric  Correction  to  Lj  RMS  (cm)  vs  Date, 
GPS  SV-9 ,  NSWC  Dahlgren,  Days  182-197  1980 
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10. 


20-MIN  IONOSPHERIC  CORRECTION  TO  E,  RMS  (cm)  VS  ELEVATION  (°) 

The  same  residuals  from  plots  7  and  8  are  plotted  versus  elevation  angle  in 
plot  10.  Here,  it  becomes  apparent  that  the  scatter  is  due  primarily  to  data 
from  low-elevation  angles.  The  SV-6  is  particularly  affected,  perhaps  because 
one  of  its  two  passes  does  not  exceed  30°. 
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Figure  A-37.  20-Min  Ionospheric  Correction  to  Lj  RMS  (cm)  vs  F.levation  (°), 
GPS  SV-5 ,  NSWC  Dahlgren,  Days  182-197  1980 
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11.  IONOSPHERIC  CORRECTION  L.  ( cm )  VS  ELEVATION  (°) 

The  ionospheric  correction  to  Lj  computed  from  two-frequency  Doppler  is 
presented  in  plot  II.  The  corrections  are  negative  during  the  portion  of  the 
pass  that  the  range  is  decreasing  and  positive  as  the  range  increases.  The 
convergence  of  the  corrections  at  the  peak  elevation  angle  is  quite  apparent  in 
the  data;  however,  it  does  not  appear  to  be  symmetric  about  zero.  Further 
investigations  are  warranted  to  explain  this  result. 
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Figure  A-41.  Ionospheric  Correction  to  Lj  RMS  (cm)  vs  Elevation  (°), 
CPS  SV-b,  NSWC  Dahlgren,  Days  182-197  1980 
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CLOCK  OFFSET,  CPS-NGRS  (nsec)  VS  DATE 

The  final  plot  demonstrates  again  the  consistency  of  the  derived  offset 
after  the  broadcast  satellite  clock  corrections  are  applied  to  the 
pseudorange  observations. 


Figure  A-45.  Clock  Offset,  GPS-NGR3  (nsec)  vs  Date 
GPS  SV-S,  NSWC  Dahlgren,  Days  182-197  1980 


Figure  A-46.  Clock  Offset,  GPS-NGRS  (nsec)  vs  Date 
GPS  SV-6,  NSWC  Dahlgren,  Days  182-197  1980 
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Figure  A-48.  Clock  Offset,  GPS-NGRS  (nsec)  vs  Date 
GPS  SV-9 ,  NSWC  Dahlgren ,  Days  182-197  1980 
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